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1 1-D Wave Equation : Physical derivation 

Reference: Guenther & Lee §1.2, Myint-U & Debnath §2.1-2.4 
[Oct. 3, 2006] 

We consider a string of length I with ends fixed, and rest state coinciding with 
X-axis. The string is plucked into oscillation. Let u {x, t) be the position of the string 
at time t. 

Assumptions: 

1. Small oscillations, i.e. the displacement u {x, t) is small compared to the length 
I. 

(a) Points move vertically. In general, we don't know that points on the string 
move vertically. By assuming the oscillations are small, we assume the 
points move vertically. 

(b) Slope of tangent to the string is small everywhere, i.e. \ux {x,t)\ <^ 1, so 
stretching of the string is negligible 

(c) arc length a (t) = -^/l + u^dx ~ 

2. String is perfectly flexible (it bends). This implies the tension is in the tan- 
gent direction and the horizontal tension is constant, or else there would be a 
preferred direction of motion for the string. 

Consider an clement of the string between x and x + Aa;. Let T (x, t) be tension 
and 9 (x, t) be the angle wrt the horizontal a;-axis. Note that 

tan 9 {x, t) = slope of tangent at {x, t) in ux-plane — — {x, t) . (1) 
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Newton's Second Law {F — ma) states that 

F - (pAx) — (2) 

where p is the hnear density of the string {ML~^) and Ax is the length of the segment. 
The force comes from the tension in the string only - we ignore any external forces 
such as gravity. The horizontal tension is constant, and hence it is the vertical tension 
that moves the string vertically (obvious). 

Balancing the forces in the horizontal direction gives 

T{x + Ax, t) cos 0{x + Ax, t)^T {x, t) cos e{x,t)^T ^ const (3) 

where r is the constant horizontal tension. Balancing the forces in the vertical direc- 
tion yields 

F = T{x + Ax, t) sine {x + Ax, t) -T{x,t) sine {x,t) 

= T {x + Ax, t) cos e {x + Ax, t) tan e {x + Ax, t) —T [x, t) cos e [x, t) tan 6* (x, t) 

Substituting (3) and (1) yields 

F = T{i&ne{x^ Ax,t) -ia.ne{x,t)) 

Substituting F from (2) into Eq. (4) and dividing by Ax gives 

dh^ _ Jt{x + Ax,t)-%{x,t) 
Pdt^^^'^>~^ Ax 

for e [x,x + Ax\. Letting Ax — > gives the 1-D Wave Equation 

= c2 = ->0. (5) 



Note that c has units fcl 



di? dx^ ' p 

1/2 



Force 
Density 



LT ^ of speed. 



1.1 Boundary conditions 



Ref: Guenther & Lee §4.2 (p. 94), Myint-U & Debnath §4.4 

In order to guarantee that Eq. (5) has a unique solution, we need initial and 
boundary conditions on the displacement u {x, t). There are now 2 initial conditions 
and 2 boundary conditions. 
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E.g. The string is fixed at both ends, u (0, t) — u{l,t) — 0, t > (homogeneous 
type I BCs) 

E.g. The string is connected to frictionless cyhnders of mass m that move vertically 
on tracks at x — 0,1. Performing a force balance at either x = or x = 1 gives 

Tsm9 = mg (6) 

In other words, the vertical tension in the string balances the mass of the cylinder. 
But r = Tcos9 = const and tan 6* — u^, so that (6) becomes 

TUx — T cos 9 tan 9 — mg 

Rearranging yields 

mg 

Ux = — , X = 0,1 

T 

These are Type 11 BCs. If the string is really tight and the cylinders are very light, 
then mg/r ^ 1 and we approximate m-j; ~ at a; = 0, 1, and the BCs become Type 
11 homogeneous BCs. 



1.2 Initial conditions 

Ref: Guenther & Lee §4.2 (p. 94) 

We need to specify both the initial position of the string and the initial velocity: 
u{x,0) = f {x) and ut{x,0) = g (x), < x < I. The idea is the same as finding 
the trajectory of a falling body; we need to know both the initial position and initial 
velocity of the particle to compute its trajectory. We can also see this mathematically. 

The Taylor series of u{x, t) about i = is 

u {x, t)=u (x, 0) + ut {x, 0) t + uu {x, 0) — + ^) 3! 

Prom the initial conditions, u{x,0) = f {x), Ut{x,0) = g {x) and the PDE gives 

Uu {X, 0) = C^Uxx {x, 0) = c^f" (x) , 

— [x, 0) = c^utxx {x, 0) = c^g" (x) . 

Higher order terms can be found similarly. Therefore, the two initial conditions for 
u {x, 0) and Ut {x, 0) are sufficient to determine u {x, t) near t = 0. 

To summarize, the dimensional basic 1-D Wave Problem with Type I BCs (fixed 
ends) is 

PDE : Uu = c^Uxx, 0<x<l (7) 
BC : u{0,t) = = u{l,t), t>0, (8) 
IC : u{x,0) = f {x) , ut{x,0) = g{x) , 0<x <l (9) 
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1.3 Non-dimensionalization 

We now scale the basic 1-D Wave Problem. The characteristic quantities are length 
L^, and time T^. Common sense suggests choosing = /, the length of the string. 
We introducing the non-dimensional variables 

X T t u{x,t) ^ f{x) T^g{x) 

J—J^ J. ^ J—J'j^ J—J^ J—J^ 

From the chain rule, 

du du dx du du du di L* du 
dx^ *d£dx^d£' di ^ *~didi^%~di 

and similarly for higher derivatives. Substituting the dimensionless variables into 1-D 
Wave Equation (7) gives 

This suggests choosing = L^/c — l/c, so that 

Uti = Uxx, < £ < 1, t > 0. (10) 

The BCs (8) become 

u (0, £) = = u (1, i) , i>0. (11) 

The ICs (9) become 

u{x,0)^f{x), Ui{x,0)^g{x), 0<x<l. (12) 

1.3.1 Dimensionless 1-D Wave Problem with fixed ends 

Dropping hats, the dimensionless 1-D Wave Problem is, from (10) - (12), 



PDE 
BC 
IC 



utt = Uxx, < X < 1 (13) 
u{0,t) ^0 ^u{l,t) , t>0, (14) 
u{x,0)^f(x), ut(x,0)^g(x), 0<x<l (15) 



2 Separation of variables solution 

Ref: Guenther & Lee §4.2, Myint-U & Debnath §6.2, and §7.1 - 7.3 

Substituting u (x, t) ^ X (x) T (t) into the PDE (13) and dividing by X (x) T (t) 
gives 

T"(t)_X"(x)_ 
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where A is a constant. The negative sign is for convention. The BCs (14) become 

u{0,t) = X(0)T(t) = 
u{l,t) = X{l)T{t) = 

which imphes 

X(0)=X(1) = (17) 

Thus, the problem for X (x) is the same Sturm- Liouville Boundary Value Problem as 
for the Heat Equation, 

X" (x) + XX (x) = 0; X{0)=X (1) = 0. (18) 

Recall that the eigenvalues and eigenfunctions of (18) are 

A„ = (nvr)^ , Xn{x) = bn sin {mrx) , n= 1,2,3,... 
The function T {t) satisfies 

T" + AT = 

and hence each eigenvalue A„ corresponds to a solution T„ (t) 

Tn {t) = an COS (nnt) + Pn sin {nnt) . 

Thus, a solution to the PDE and BCs is 

Un (x, t) = (an COS (mrt) + /3n sin (mrt)) sin {mix) 

where we have absorbed the constant 6„ into a„, 

In general, the individual {x, t)'s will not satisfy the ICs. Thus we sum infinitely 
many of them, using the principle of superposition. 



u {x, t) = Un {x, t) = {an COS {nnt) + /?„ sin {mrt)) sin {mrx) 

n=l n=l 

Imposing the ICs gives 

oo 

f {x) — u{x,0) — an sin {mrx) 

n=l 

oo 

g {x) — Ut{x,0) — nnPn sin {nnx) 

n=l 

Therefore, the coefficients q;„ and /?„ are given by 

an — 2 / f {x) sin {mrx) dx 



2 

Pn = — / 9 (x) sin {mrx) dx 
nn Jo 
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Note: The convergence of this series is harder to show, because we don't have 
decaying exponentials e~"^'^^* in the sum terms (more later) . 

Note: Given BCs and an IC, the wave equation has a unique solution (Myint-U 
& Debnath §6.3). 



3 Interpretation - Normal modes of vibration 

[Oct 5, 2006] 

Ref: Guenther & Lee p. 100 problem 5 
The terms 

Un {x, t) = {an COS {niTt) + sin {mrt)) sin {uttx) 

for n — 1, 2, 3, ... are called the normal modes of vibration . The solution u {x, t) is a 
superposition of the normal modes Un {x, t). In physical variables, the normal modes 
are 

, , , / /nvrc A . fmic ,\\ . f mix' 

Un[x ,t) ^ [an cos y-j-t j + Pn sm y-^t ))^^^\ ~1~ 

3.1 Frequency and Period 

A function / (t) is periodic if for some real number T, 

f(t + T)^f(t) 

for all t. If t measures time (either physical or dimensionless) , we define the period 
of a function as the smallest number T such that f (t + T) = f (t). We say that / {t) 
has period T or, equivalently, that / (t) is T-periodic. The period T has the same 
dimensions as t. If t is dimensionless then so is T; if t has the dimensions of time, 
so does T. Note that each normal mode Un{x,t) has period 2/n, and in physical 
variables, w't i^'^t') has period 21/ (nc). 

The analog to period for spatial coordinates is the wavelength. The wavelength of 
a function g (x) {x is a scaled or physical spatial coordinate) is defined as the smallest 
L such that g {x -\- L) = g{x). We say g (x) is L-periodic. Again, the wavelength 
L has the same dimensions as x. If x is dimensionless then so is L; if x has the 
dimensions of length, so does L. 

The frequency / is defined as / = 1/T = 1/period. In physical variables, each 
normal mode u'n {x', t') has frequency 

_ujn _nc _ n Ft 
~ 2n ~ Yl ~ 21V ~p 
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with dimensions 1/time (e.g. Hz = cycle/sec = sec~^). 

The angular frequency is defined as a; = 27r/ where / is the frequency. Each mode 
has angular frequency a;„ = 27r/„ = mrc/l. In terms of the frequency, we can write 
the normal mode as 



U^{x,t) = [anCOSy— j + PnSm\^—i-t j j SI 

= {an COS {2nfnt') + (5n sin (27r/„t') ) sin 
— {an cos {uJnt') + Pn sin {uJnt')) sin 



sm 



nvrx 



mix 



mrx 



I 

The first harmonic is the normal mode of lowest frequency, ui {x, t) or in physical 
variables, u[ {x',t'). 

The fundamental frequency is /i, i.e. the frequency of the first harmonic. In 
dimensionless variables, /i = 1/2. In physical variables, f[ — cj {21) — (^\/t/pJ /2l. 

Note that = n/i, in other words, the frequencies of higher harmonics are just 
integer multiples of the fundamental frequency /i. 

Examples. Check for yourself that sin (ttx) is 2-periodic i.e. has period 2, sin {x) 
is 27r-periodic, sin(n7ri) is 2/n-periodic, sva.{m:ct/l) has period 21/ {nc), sm{mrx/l) 
has period 21 /n. Note that if a function has period T, then f {t + mT) — f {t) for all 
m = 1, 2, 3, .... In particular, 

sin {mrx) , cos {mrx) , sin (nvrt) , cos {mrt) 

are all 2/n-periodic, and hence 

2m\\ 

sin \ nn \ X -\ ] ) = sin mrx, cos \nn \x-\ j J ^ ^'^'"'^ 

/ / 2m\\ 
smnnt, cos I nvr I t H I I = cosnnt 

ire, the dimensionless solution u {x, t) of the wave 
equation has time period 2 {u {x,t + 2) = u {x, t)) since 

OD OO 

u {x, t) = Un {x, t) = {an COS {mit) + (3n sin (nvrt)) sin {mix) 

n=l n=l 

and for each normal mode, Un{x,t) = Un{x,t + 2) (check for yourself). Thus, the 
period of u{x,t) is the same as that for the first harmonic Ui{x,t). In physical 
variables, the period of u' {x\t') is T[ = 21/ c = 21/ \Jt / p and the frequency is = 
\/T[ = c/{2l). 





2m 


[x + 






n 




2m 


(■* 






n 


1,2, 


3... 
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3.2 A note on dimensions 

I've tried to denote dimensional quantities with primes or specifically comment whether 
or not we are working with physical x, t (a length and time) or with dimensionless x, 
t (dimensions 1). However, in general you should always ask the question, "What are 
the dimensions?" The quantity c/l has dimensions of 1/time since I have defined c, / 
to be a speed and a length. The argument of any mathematical function like cos, sin, 
exp, etc. must be dimensionless. The cosine of 1 m or 2 s does not make sense. Thus, 
the quantity t in cos [mit) is dimensionless. Note that radians are a dimensionless 
measure of angle. Consider the quantity t in cos {uit) or cos (27r/t) or cos (^) , where 
you are given that uj and / have dimensions of 1/time, c is a speed, and I is a length. 
Then t must have dimensions of time, so that the arguments ujt, 27r ft, or irct/l of 
cosine are dimensionless. 

3.3 Amplitude 

Note that ti„ {x, t) can be written 

Un {x, t) — 7„ sin (mrx) sin {mrt + ipn) (19) 

where 7„ = ^/Q^~\~|3f^, '0n — arctan (^^^- At each point x, the n'th mode vibrates 
sinusoidally according to 

Un {x, t) = An (x) sin {mrt + ipn) 

where 

An {x) = 7„ sin {wKx) . 

The mode it„ (x, t) vibrates sinusoidally in time between the two limits ±A„ [x). We 
call A {x) the time amplitude of the mode Un {x,t), since this sets the bounds on the 
oscillations in time. Locations where An (x) — are called nodes and locations where 
\An {x)\ — 7„ are called antinodes. 

4 Conservative system and energy 

Ref: Guenther & Lee p. 102, Myint-U & Debnath §4.12 (problem 28 p. 97, 1987) 

Recall the solution to the Heat Problem with a homogeneous PDE (i.e. no sources, 
sinks) Homogenous Type I BCs {u — Bit x — 0,1) is 

oo 

u {x, t) = ^2 sin (nnx) e"'*^'^^* 

n=l 
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and as i — > oo, u {x, t) — > 0. Thus, the rod looses heat energy and with it the memory 
of the initial state, or initial condition. 

In contrast, the solution to the wave equation with Homogeneous Type I BCs 
(fixed ends, u — a,t x — 0,1), 

oo 

u {x, t) — {an COS (mrt) + /3„ sin (mrt)) sin (mrx) 

n=l 

The oscillations do not decay, since we ignored gravity and resistive forces. Once 
plucked, the string vibrates/oscillates forever, with period 21 /c in physical coordi- 
nates. We call this type of system conservative: energy is conserved. In addition, the 
system comes back to its initial condition periodically - i.e. it maintains a memory 
of its initial state. 

We define the energy of the system as 

Total energy E' {t') = PE' {f) + KE' (t') 

where primes denote physical or dimensional variables. The local kinetic energy is 

KE'{x\t')^\mv'^\{pl^x) 

and the local potential energy is 

PE' {x' ,t') = work done getting to displacement u 
— Force x change in length of string 

Making a displacement Au of a segment of string from x to x + Ax results in a change 
in the string segment length, initially of length Ax, of 



fAuV 



Al^ \/{Axf + {Auf -Ax^Ax { -1 



Recall the binomial expansion ^/T+~a? = 1 + ^0^ + (a'^). Then 



1 f AuV . ^ f f AuVW AxfAu^"^ 



Ai = Aa; l+- — -1 + 0,,, ,,, 

2 \Ax) WAxJ I I 2 \Ax 

and hence the potential energy is 
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as Ax — > 0. Thus, the total energy is, in dimensional form, 

where we use the shorthand notation u^, = {du/dt'Y . Substituting (? — t j p and the 
dimensionless variables into (20) gives the dimensionless total energy 

EiS) = ^ = \\\u\^ul)dx (21) 



rl 2 Jo 

4.1 Conservation of energy for a single mode 

[Oct 17, 2006] 

The energy of a given normal mode n„ (x, i) is found by substituting the form (19) 
into Eq. (21) for the energy, 

2 Jo 



cos^ (nvrt + ipn) J sin^ (nvra;) dx 
+ sin^ (nTrt + ipn) J cos^ (nvra;) (ix 



Note that 



sin [nTTX) = 7: ~ 7^ cos {2mrx) , cos [mix) = 7^ + 7: cos {2n7rx) 

so that 

^ 1 /"^ 1 

sin^ (riTTx) dx — / cos^ (nvrx) dx — - 

2' 7o 2 

Substituting these integrals into the energy (t) gives 

£;„ (i) = i!^!!^ (cos^ (riTTi + V'n) + sin^ (nTri + V^„)) = 
Thus the total energy (kinetic + potential across the rod) of each mode is constant! 

4.2 Energy conservation for multi-mode wave 

To show the energy for a multi-mode wave is conserved, we proceed as follows. Dif- 
ferentiating the dimensionless energy equation (21) gives 

10 



Replacing Uu — u^x (the wave PDE) gives 

dE _ 
~dt 



: 11,11 rr + UxUxt) dx 
{utUx)^ dx 

= K«Ji=o (22) 
Differentiating the BCs in time t gives 

^«(0,t)=0, ju{l,t) = Q 

Therefore, (22) becomes 



Thus 



E{t) = const = E{0) = ^J^' {(/' {x)f + {g {x)f] dx. 



Thus the energy of a multi-mode wave system remains constant in time. This is 
not surprising, since all the forces involved are conservative; we have not included 
damping or dissipative forces. 

An alternate method to show E (t) — const employs the infinite series solution 
derived using the method of separation of variables: 

2^2^ dt dt dx dx I 

^n=l m=l / 
^ f dUn dUm _^ dUn dUr„ . 
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2 ^ ^ /n V dt dt dx dx 

n=l m=l 



Substituting for the normal modes ti„ {x, t) = 7„ sin {mrx) sin {mrt + ipn) gives 
E (t) — - jnlmnmn'^ cos {nnt + ipn) cos {mnt + ipm) 

n=l m=l 

X / sin in-Kx) sin {rxi'nx) dx 
Jo 

+- ^ ^ ^nlmnmn'^ sin {nnt + ipn) sin {mnt + ipm) 

n=l m=l 

/ COS (nnx) cos (mnx) dx 
Jo 



n=l m=l 

X 
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Substituting the orthogonality relations gives 

^ ^ XI XI Inlmnrm:'^ cos {n-Kt + V'n) cos (mTrt + V'm) ^ 

n=l m=l 

+^ X X 'ynlmnmn'^ sin (nvri + sin (rnvri + -0^) ^ 

n=l m=l 

= I X ("^n^^)^ (^°^^ (^^^ + '^") + ^^^^ (^^^ + '^ri)) 
n=l 

?i=l 

which is constant for all t. 

5 D'Alembert's Solution for the wave equation 

Ref: Guenther & Lee §4.1, Myint-U & Debnath §4.3 

5.1 Motivation 

Note that each normal mode can be written in alternative form: 

Un {x, t) — {an COS (mrt) + /3„ sin (mrt)) sin (mrx) 

— - (an sin (riTT {x — t)) — cos {nn {x — t))) 

1 

+ - (a„ sin (riTT {x + t)) + cos (wtt {x + t))) 

where we used the trig identities 

sin (a + 6) + sin (a — &) = 2 cos h sin a 
cos (a + 6) — cos (a — &) = 2 sin a sin 6 

Therefore, the mode Un [x-, t) is the sum of a function of {x — t) and a function of 
(x + 1). In physical variables, we can write the solution as the sum of a function of 
x' — ct' and a function of x' + ct'. 

5.2 Change of variable 

We can define the new coordinates 

^^x-t, r]^x + t (24) 
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and let v rf) — u (x, t). The chain rule implies that 

du dv dv dr] dv dv 

dx dx dr] dx dr] 

du dv d^ dv dr] dv dv 
dt d^ dt dr] dt d^ dr] 



Similarly, 



d'^u d dv d dv / d'^v d'^v \ / d'^v d'^v 
dx^ dx d^ dx dr] \d^'^ d^drj J \d^dr] dr]'^ 



d'^v d^v d'^v 

= h2 \ 

d^dr] dr]"^ 

d'^u / d'^v d'^v \ f d'^v d'^v 

dt"^ \ d^'^ d^dr] ) \ d^dr] dr]"^ 

d'^v d'^v d'^v 

= 2 \ 

d^dr) dr)'^ 

Substituting these derivatives into the PDE uu = Uxx yields 

d'^v d'^v d'^v d'^v d'^v d'^v 

h 2 \ = 2 \ 

d^'^ d^drj drj'^ d^'^ d^dr] dr}'^ 

Simplifying and dividing by 4 gives a new form of the wave equation, 

d'^v (C, v) 



(25) 



d^drj 

5.3 Forward and backward waves 

We can write the new form (25) of the wave equation in two ways: 

dfdv{^,r])\ ^ d f dv{^,r]) \ 



d^ \ dr] J ' dr] \ d^ 
Integrating the first equation in ^ gives 

dv (^, 77) 



dr] 



G [r]) (26) 



for some function G {r]) (due to partial integration with respect to ^). Define the 
antiderivative of G [r]) as Q (77), i.e. such that Q' {vj) = G {vj). Substituting this into 
(26) and rearranging gives 
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Integrating the right equation in rj yields 

where the function P {^) is due to the partial integration with respect to rj. Therefore, 
we can write 

Substituting for r] and ^ from (24) and recalling that u {x, t) = v rj) gives 

u{x,t)=v{C,r]) = P{x-t) + Q{x + t) (27) 

5.3.1 Forward wave 

The forward wave is the function P {x — t) which represents a wave travelling in 
the positive x-direction with scaled velocity 1. In physical coordinates, the function 
depends on x — ct and the speed of the wave is c = \Jt j p. The shape of the wave is 
determined by the function P [x) and the motion is governed by the line x — t = const. 
The wave (same shape) moves forward in time along the string. 

How do we find the speed from the fine x' — ct' — h — const! Well, in time At', 
we've moved 

Ax' ^{h + c {t' + At')) - ih + ct') = cAt' 
The speed is distance traveled over elapsed time: Ax'/ At' = c. 

5.3.2 Backward wave 

The backward wave is the function Q {x + t) which represents a wave travelling in the 
negative x-direction with scaled speed 1. The wave shape is determined by Q (x) and 
the value of the wave is constant along the lines x + t — const (in physical variables, 
x' + ct' — const and speed is c) . 

5.4 Characteristics 

Ref: Myint-U & Debnath §3.2(A) 

The solution to the wave equation is the superposition of a forward wave P {x — t) 
and a backward wave Q {x + t), both with speed c. The lines x±t — const are called 
characteristics. 
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5.5 Determining the shape functions 

The shapes of the forward and backward waves, P [x) and Q (a;), are determined from 
the initial conditions, 

u{x,Q) = f{x)^P{x) + Q{x), (28) 
Ut{x,Q) = g{x)^-P'{x) + Q'{x) (29) 

To obtain the second equation, the chain rule was used: 

du 
di 
' d 



ut (x, 0) 



{x,t) 



dt 



{P{x-t) + Q{x + t)) 



t=0 



= [-P' {x-t) + Q' {x + 1)],^, 
= -P'(x) + Q'(x) 

It is important to differentiate in time first and then set t = (that is what Ut {x, 0) 
means!). Now, integrating (29) in x from to x gives 

Q{x)-P {x) - {Q (0) - P (0)) = r g is) ds (30) 

Eqs. (28) and (30) can be solved for P (x) and Q (x): 

Qix) = ^(^f{x) + £g{s)ds + QiO)-PiO) 



fix)- / g{s)ds-Q{0) + PiO) 



P{x) 



5.6 D'Alembert's solution to the wave equation 



(31) 
(32) 



[Oct 19, 2006] 

Ref: Guenther & Lee §4.1, Myint-U & Debnath §4.3 

Summarizing our results: from (27), (31) and (32), the wave equation and initial 
conditions 



u {x, 0) 
Ut {x, 0) 



fix) 

g{x) 



has the solution 

u(x,t)^P(x-t) + Q(x + t)^l[f(x-t) + f(x + t)] + l-[^ g(s)ds (33) 

This is called D'Alembert's solution to the wave equation. We have not yet considered 
the BCs [we will soon], so right now we're thinking about an infinite string. 
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5.6.1 Solution method and domain of dependence 

At a given position x = xq on the string, the solution at time t = to is 

rxo+to 

u {xq, to) = - [/ {xo -to) + f {xo + to)] + - g{s) ds. 

In other words, the solution is found by tracing backwards in time along the charac- 
teristics X — t — Xq — to and x -\- 1 — Xo -\- to io the initial state (/ (x), g (x)), then 
applying (33) to compute u (xq, ^o) from the initial state. Information from the initial 
state from the interval xo — to < x < xo -\- to is, all that is needed to find u (xq, to). In 
the te-plane, we can think of a triangle opening backwards in time from {xo,to) to 
the line t — Q. Compare this to the rectangle {0 < x < 1, < i < T} we required to 
find the solution u (a:, T) for the heat equation. 

For simpler forms of / {x) and g [x) the solution method can be simplified signif- 
icantly, and only a handful of characteristics will be needed. 

5.6.2 Simplified solution method for infinite string 

For the examples we'll be considering, / [x) and g {x) will be case functions. The 
general approach is then: 

Step 1. Write down D'Alembert's solution, 

= + + + ^ 9{s)ds 

Step 2. Identify the regions. In general, the function / {x) and g {x) are case 
functions. You need to determine various regions by plotting the salient characteristics 
x ±t = const. The regions determine where x — t and x + t are relative to the cases 
for the functions / {x) and g {x) and tells us what part of the case functions should 
be used in each region. 

Step 3. Determine the solution in each region. 

Step 4. For each specific time t = to, write the a;-intervals corresponding to the 
regions. 

For example, consider / (x) and g (x) of the following form, 

, iF(x), \x\<l {g(x), \x\<l 

I 0, |x| > 1 I 0, \x\ > 1 

Step 1. D'Alembert's solution to the wave equation is 

"X+t 



^i^,t) ^^[f {x -t) + f {x + t)] + ^ ^ g{s)ds 
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Figure 1: Regions of interest separated by four characteristics. 

Step 2. Identify the regions. The functions / (x) and g ix) arc equal to functions 
F [x) and G[x\ respectively, for < 1 and are zero for |x| > 1. Thus, the regions 
of interest arc found by plotting the four characteristics a; ± t = ±1 (Figure 1). The 
regions are identified in the plot, and are given mathematically by 





{{x,t) 




l<x-t<l 


and -l<x + t<l} 




{{x,t) 




l<x-t<l 


and x + t > 1} 


Rz = 


{ix,t) 


X 


— t < —1 and 


-l<x + t<l} 


R4 = 


{{x,t) 


X 


— t < — 1 and 


x + t>l} 


R5 = 


{(x,t) 


X 


+ t<-l}, 




Rq — 


{{x,t) 


X 


-t>l} 





The regions determine where x — t and x + t are relative to ±1, which tells us what 
part of the case functions / {x) and g {x) should be used. It is helpful to define the 
fines 

XA{t) ^ -t-1, XB{t)^t-l, Xc{t)^-t + l, XD{t)^t+l. 

Step 3. Consider the solution in each region. In Ri, we have |a; ± t| < 1, so that 
(34) implies 

rx+t rx+t 

f{x-t)^F{x-t), f{x + t)^F{x + t), / g{s)ds^ G (s) ds, 

J x—t J x—t 
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and hence 



1 1 r^"^ 

u(x,t) ^ - (F (x -t) + F (x + t)) + - G (s) ds. 

^ ^ Jx-t 

In region i?2, we have — 1 < a; — t < 1 and a; + t > 1 so that 



/(x + i) = 0, f{x-t)^F{x-t), 



rx+t 

/ 9 (s) ds 

Jx-t 



9 (s) ds 



x—t 



x+t 



9 {s) ds 



G (s) ds + 0, 



x—t 



and hence 



, , F{x-t) 1 , , ^ 



x—t 



In region R^, we have —l<x + t<l and x — t < —1 so that 



/(a: + t) = F(a; + i), / (x - = 0, 
and hence 



x+t 



9 {s) ds 



x—t 



x+t 



G (s) ds, 



u {x, t) 



F{x + t) 1 



+ 



x+t 



2 2 

In region R^, x — t < —1 and x + 1 > 1, so that 

f{x + t)^O^f{x-t), 



G (s) ds. 



rx+t pi pi rx+t pi 

/ g{s)ds= I g{s)ds+ I g{s)ds+ l g (s) ds = + G (s) ds + 0, 

Jx-t Jx-t J-1 Jl J-1 



G (s) ds — const. 



and hence 

1 r+* 1 

^) = 2 J ^ 9{s)ds^- 

In regions R5 and i?6, / (a; + t) = = / (x — t) and (s) = for s e [x — t,x + t], 
hence u — 0. To summarize, 

' ^{F{x-t) + F{x + t)) + ^ J^^^G{s) ds, {x,t) e Ri 

^ + Ux-tGis)ds, {x,t)eR2 

u{x,t) = \ ^ + iJ_-+*G'(s)ds, ix,t)eRs (36) 

I /^^ G (s) (is, (x, t) e i?4 

{x,t)ER5,R6 

Step 4. For each specific time t — to, write the x-intervals corresponding to the 
intersection of the sets i?„ with the line t — to. This amounts to computing the values 

of xa {t), xb {t), Xc (t), xd {t) for each time 



t 





1/2 


1 


2 


XA (t) = -l-t 


-1 


-1.5 


-2 


-3 


XB{t)=t-l 


-1 


-0.5 





1 


Xc (t) = -t + 1 


1 


0.5 





-1 


xnit) =t + l 


1 


1.5 


2 


3 



(37) 
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At t — 0, we use Table (37) and Figure 1 to find the x intervals i?^ corresponding 
to the intersection of i?„ with the line t — 0: 

R', = {-oo,-l], R[ = [-l,l], R'^^[l,oo). 

In Ri, we have (recall that t — 0), 

1 1 f^~^^ 

u(x,0) ^ - (F (x - 0) + F (x + 0)) + - G (s) ds ^ F (x) 

2 2 J^_Q 

Similarly, we can check that in the other regions, u — 0, so that (36) becomes 



F(x), XG /?; = [-!,!] _ 1 




a; 


0, xeR'r^UR'e 1 


I 0' 


a; 



At t = 1/2, we use Table (37) and Figure 1 to find the x intervals i?^ corresponding 
to the intersection of i?„ with the line t — 1/2: 

i?5 = (-00,-1.5], i?i = [-0.5,0.5] , i?6 = [l-5,oo), (38) 

i?3 = [-1.5, -0.5] , i?2 = [0.5, 1.5] . (39) 

At t = 1, we use Table (37) and Figure 1 to find the x intervals i?^ corresponding 
to the intersection of i?„ with the line t = 1: 

i?5 = (-00,-2], i?3 = [-2,0], i?2 = [0,2], i?6 = [2,oo). (40) 

At t = 2, we use Table (37) and Figure 1 to find the x intervals i?^ corresponding 
to the intersection of Rn with the line t = 2: 

i?5 = (-00,-3], i?3=[-3,-l], i?4=[-l,l], i?2 = [l,3], i?6 = [3,oo), 

(41) 

At this point, you usually want to draw the wave profiles for certain times t — to 
(see problem set 3 for more). 

5.6.3 Worked example 

For an infinitely long string, consider giving the string zero initial displacement 
u {x, 0) = and initial velocity Ut {x, 0) = g (x). Suppose 



g{x) 



cos^ (fx) , -1 < X < 1 
otherwise 



The ICs have the form considered above for F (x) = and G (x) = cos^ (f^)- 
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Step 1. D'Alembert's solution (33) becomes 



U (x, t) 



g (s) ds 



x—t 



Step 2. The regions are the same as those in (35) and are plotted in Figure 1 
above. 

Step 3. Determine u {x,t) in each region. From (36), we have 



u {x, t) 



' \^l-^lG{s)ds, {x,t)eR, 

'^Jl^G{s)ds, {x,t)eR2 

lj:^'G{s)ds, {x,t)eR3 

\ \\ G (s) ds, (x, t) e i?4 

{x,t)eR5,R6 



Note that 



dx 



(1 + COS (nx)) dx 



b — a sin {nb) — sin (vra) 



+ 



27r 



Thus 



rx+t 

Jx-t 
1 



^, , , sin (n (x + 1)) — sin (n (x — t)) 1 / s . / s 

^(5)^5 = i H ^ = i + - cos (nx) sin (tt^) 



G (s) ds 



x—t 
x+t 



G (s) ds 



27r 

1 — (x — t) sin (tt (x — t)) 

2 2^1^ 
1 + a; + t ^ sin (tt (a; + t)) 



27r 



G{s)ds = 1 



TT 



Thus 



u (x, i) = < 



I + ^ sin (tt^) cos (ttx) , (x, i) e i?i 

1— (a;— t) sin(7r(a;— t)) 

4 47r ' 

1+x+t I sin(7r(j;+t)) 



4 + 47r 
1 

2' 





(x, t) e i?2 
(x, i) e i?3 
(x, i) e it:4 
{x,t) e R^.Rq 



(42) 



Step 4. We consider early, intermediate and later times, t — 1/2, 1, 2. Ati = 1/2, 
the regions Rn are given by (38), (39) and (42) becomes 



u X 



X 



cos TTX 

cos nx 



+ 4. 



4-n- ' 
cosTra; 



l<x<l 

-i < x < 



1 



— 2 

ki>i 
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Figure 2: The profiles of tlie displacement u{x, t) for the times t = 0, 1/2, 1, 2. 



At t — 1, the regions Rn are given by (40) and (42) becomes 



u {x, 1) 



^ + ^, 0<x<2 

2^_sjn^^ -2<X<0 







Ixl > 2 



At i = 2, the regions Rn are given by (41) and (42) becomes 

r 3^_sH^ i<x<3 

4 47r ' — — 

3+x 1 sm(7ra;) _ 

M (X, 2) = < 4 4^ 



1 
2' 





, -3<x<-l 
-1 < a: < 1 
\x\ > 3 



In Figure 2, the profiles of the displacement u{x,t) are plotted for the times t 
0,1/2,1,2. 



6 Waves on a finite string 

[Oct 24, 2006] 

Ref: Myint-U & Debnath §4.4 - 4.6, Guenther & Lee §4.5 

We now consider D'Alembert's solution for a finite string. The dimensionless 
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problem with homogeneous Type I BCs (ends of string fixed) is 

Utt = Uxx 

u{0,t) = Q^u{l,t) (43) 

u{x,0) = fix) 

ut{x,0) = g{x) 

We found that this has solution 

oo oo 

u {x, t) = Un {x, t) = {an COS {nnt) + /3„ sin {nnt)) sin {nnx) 

n=l n=l 

where 

an — 2 / f (x) sin (nnx) dx 
Jo 

2 

(5n = — / g ix) sin (nTTx) dx 
mr Jo 

We wrote, equivalently, that 

u{x,t) ^ P {x - t) + Q {x + 1) 

where 

^ oo 

P {x — t) — - («„ sin {mr {x — t)) — /3„ cos {mr {x — t))) 

n=l 
^ oo 

Q{x + t) = -^(Q;„sin(?i7r(a; + i)) +/3„cos(n7r(a; + t))) 



2 

n=l 



6.1 Zero initial velocity 

The simplest ICs involve the case of zero initial velocity i^t (x, 0) =0 and a specified 
initial displacement u (x, 0) — f (x). In this case, /3„ = for all n and 

oo 

u {x, t) = an COS {mrt) sin (riTrx) = P (x — i) + Q (a; + i) 

n=l 

where 



1 1 - 

P (s) = Q (s) = - J] an sin riTrs = -/ (s) (44) 

n=l 

where /(x) is the odd 2-periodic extension (Sine series) of / (x). We can thus write 

u{x,t) = P{x-t) + Q{x + t) = ^ (/ {x-t) + f {x + t)) . (45) 
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This looks like D'Alembert's solution for the infinite string - but / replaces / (x) - 
so we'll see that some properties of the solution are different. 

To see directly that (45) satisfies the BCs, note that since / is odd, / {—t) — —f (t), 
and since / is 2-periodic and odd, 

/(l-t) = /(-l-t) = -/(l + t) 

Thus, 

uiO,t)^^(^f{-t) + fit))-0 

Therefore, (45) satisfies the Type I BCs (fixed ends). One can also check directly 
that (45) satisfies the ICs and PDE. 

To summarize, we have written the displacement of the string in terms of the odd 
2-periodic extension of the initial condition / (x). Eq. (45) is the solution for the finite 
string problem (43) with Type I BCs (fixed ends), zero initial velocity Ut{x,0) — 
and initial displacement u {x, 0) — f (x). 

To draw the solution, note that x only goes from to 1, but time t G [0, oo). Thus 
X + t e [0, oo) and x — t G [0, — oo). We use the 2-periodic and odd properties of / to 
see what u{x,t) looks like for < a; < 1. 

Example. Suppose f (x) is a thin pulse symmetric about x = 1/2 and with a 
maximum at x = 1/2. The initial pulse breaks into forward and backward waves, 
each propagating with speed 1 in opposite directions. The peak of the forward wave 
(at a; = 1/2) first reaches the end of the string x = 1 when 

^ _ _ 1 1 

~ 2 ^ ~ 2 

At the same time, the peak of the backward wave reaches the end x = 0. Then what 
happens? On the first half of the string < x < 1/2 and times 1/2 < t < 1, 

-l<x-t<0 

and thus 

f(x-t)^-f{t-x)^-f{t-x) 

Recall that / (x) is the shape of the original pulse, so that on the first half of the 
string we have the negative of the forward wave. In other words, the forward wave 
went off the end x — 1 and reappeared a,t x — 0, UPSIDE DOWN (i.e. with negative 
sign)! On the second half of the string 1/2 < x < 1 and times 1/2 < i < 1, 

l<x+t<2 
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and thus 

f + t) ^ -f + t - 1) ^ -fix + t - 1) 

Thus the backward wave went off the end a; = 0att = f/2 and reappeared at the 
other end x — 1 (also ai t — 1/2) and upside down (negative sign). At t — 1, the 
waves meet in the center, 

u (x, 1) = I (/ {x-l) + f (x + (-/ (x) - f (x)) = -/ (x) = -/ (x) 

and the displacement is the upside-down version of the initial displacement. The 
process continues in this manner. The characteristics x ±t crisscross each other in 
the xt-plane. We can summarize the behavior using the following table: 



time 


X range 


X — t 
range 


x + t 
range 


left half string 


right half string 


< i < 1/2 


[0, 1/2] 




[0,1] 


bkwd wave 






[1/2, 1] 


[0,1] 






fwd wave 


1/2 < f < 1 


[0. 1/2] 


[-1.0] 




fwd wave, 

upside down 






[1/2, 1] 




[1,2] 




bkwd wave, 
upside down 
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